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Abstract 

o ; 

^ We study the g-essence model with Yukawa interactions between a scalar field and a Dirac 

field ^. For the homogeneous, isotropic and flat Friedmann-Robertson-Walker universe filled 
with the such g-essence, the exact solution of the model is found. Moreover, we reconstruct 
the corresponding scalar and fermionic potentials which describe the coupled dynamics of the 
l/^ . scalar and fermionic fields. It is shown that some particular g-essence models with Yukawa 

' interactions correspond to the usual and generalized Chaplygin gas unified models of dark 

energy and dark matter. Also we present some scalar-fermionic Dirac-Born-Infeld models 

O corresponding g-essence models with Yukawa interactions which again describe the unified 
' dark energy-dark matter system. 

o 
U 

1 Introduction 

^ ■ One of the most puzzling discovery of the last years in physics is the current acceleration of the 

0^ I universe [T]-[2]- An unknown energy component, dubbed as dark energy, is proposed to explain this 

^>0 ■ acceleration. Dark energy almost equally distributes in the universe, and its pressure is negative. 

The simplest and most theoretically appealing candidate of dark energy is the cosmological constant 
that is the ACDM model. In this case, the equation of state parameter w = — 1. Although the 
I ACDM model is in general agreement with the current astronomical observations, but has some 

difficulties e.g. to reconcile the small observational value of dark energy density with estimates from 
quantum field theories. So although the ACDM model is the most obvious choice, but it suffers 
from coincidence problem and the fine-tuning problems. It is thus natural to pursue alternative 
possibilities to explain the mystery of dark energy. In order to explain the acceleration that is 
. dark energy, many kinds of models have been proposed, such as quintessence, phantom, k-essence, 

H I tachyon, f-essence, Chaplygin gas and its generalizations, etc. 

In the last years, the k-essence model has received much attention. It was originally proposed as 
a model for inflation [3], and then as a model for dark energy [l]-[7]. Since from it was proposed, k- 
essence was been studied intensively. It is still worth investigating in a systematic way the possible 
cosmological behavior of the k-essence. Quite recently, the so-called g-essence model has been 
proposed |13| . which is a more generalized model than k-essence. In fact, the g-essence contains, 
as particular cases, two important models: k-essence and f-essence. Note that f-essence is the 
fermionic counterpart of k-essence. 

To our knowledge, in the literature there are relatively few works on the dark energy models 
with fermionic fields. However, in the recent years several approaches were made to explain the 
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accelerated expansion by choosing ferinionic fields as the gravitational sources of energy (see e.g. 
refs. |10]-|36|). In particular, it was shown that the fermionic field plays very important role in: 
i) isotropization of initially anisotropic spacetime; ii) formation of singularity free cosmological 
solutions; iii) explaining late-time acceleration. In particle physics, the Yukawa interaction plays 
an important role. It has the form 

U = (1.1) 

It describes the interaction between a scalar field and a Dirac field ip. Some properties of the 
Yukawa interaction (1.1) related with the gravitational field were considered in [1I]-[1S]. With the 
Yukawa interaction (1.1) is related the so-called Yukawa potential, 

V{x) = -^^.T-ie-™^ (1.2) 

which is negative, that is, the corresponding force is attractive. The relation between the Yukawa 
potential (1.2) and the accelerated expansion of the universe were studied by some authors (see e.g. 
[35]). In this paper, we focus on so-called g-essence model [13] which is some hybrid construction of 
k-essence and f-essence. If exactly, we will consider the g-essence with the Yukawa interaction (1.1). 
The formulation of the gravity- fermionic theory has been discussed in detail elsewhere |37)-|40j.. 
so we will only present the result here. 

This paper is organized as follows. In the following section, we briefly review g-essence. In Sec. 
3, we introduce the g-essence model with the Yukawa interaction. In Sec. 4, we construct the so- 
lution of the particular g-essence model with the scalar-fermionic Yukawa interaction. The unified 
scalar-fermionic Chaplygin gas model of dark energy and dark matter from the g-essence model 
with the Yukawa interaction were constructed in section 5 and its extension for the generalized 
Chaplygin gas case in section 6. The scalar-fermionic Dirac-Born-Infeld (DBI) counterpart of the 
g-essence model with the Yukawa interaction (3.1) was constructed in section 7. Finally, we shall 
close with a few concluding remarks in Sec. 8. The metric signature used is (-I-, — , — , — ) and units 
have been chosen so that SttG ~ c ~ h ^ 1 . 



2 Basics of g-essence 

The action of g-essence has the form [T5] 

S = [ d^Xy^[R + 2K{X,Y,(l},^,'4>)], (2.1) 



where K is the g-essence Lagrangian and is some function of its arguments, is a scalar function, 
ip = (V"!, V'2) '03j "04)"^ is a fermionic function and 'ip = ij}^"f^ is its adjoint function, the curvature 
scalar R. Here 

X = o.Ss^'^v^^v,^, y = o.5i[^r^i?^v-(^pV^)r^V'] (2.2) 

are the canonical kinetic terms for the scalar and fermionic fields, respectively. and I?^ are the 
covariant derivatives. The fermionic fields are treated here as classically commuting fields. The 
model (2.1) admits important two reductions: k-essence and f-essence. In this sense, it is the more 
general essence model and in [13J it was called g-essence. Note that to find the equations of motion 
we need the variations 

SV^ = -0.5.g,,,^/^5g'^^ (2.3) 
5R = (i?^, + .g^.D - V^V,)(5g'^^ (2.4) 
SK = K'S(j) + K"Si: + K"'Si:, (2.5) 

where Vj^y^ = di.Vf^—V^^.Va is the covariant derivative of a vector and the curved d'Alembertian 
on a scalar (jj is 

^9^(yq^a^0). (2.6) 

We now consider the dynamics of the homogeneous, isotropic and flat FRW universe filled with 
g-essence. In this case, the background line element reads 

ds^ ^ - [dx^ + dy"^ + dz^) (2.7) 



2 



and the vierbein is chosen to be 

{e'^) = diag{l,l/a,l/a,l/a), (e'^J = diag(l,a,a,a). (2.8) 
In the case of the FRW metric (2.7), the equations corresponding to the action (2.1) look like 

m 





3ij2 - p 


= 0, 


(2.9) 




2H + 3H'^ +p 


= 0, 


(2.10) 






= 0, 


(2.11) 


KyTp + 0.5{3HKy - 


h kY)i^ - ij^'K^ 


= 0, 


(2.12) 


Kyi! + 0.5{3HKy - 




= 0, 


(2.13) 




p + 3Hip + p) 


= 0, 


(2.14) 



where Kx = dK/dX, Ky = dK/dY, = dK/dcj), = dK/d%l), = dK/d%l> and 
H = aja denotes the Hubble parameter, the dot represents a differentiation with respect to time 
t. Here the kinetic terms, the energy density and the pressure take the form 

X = O.50^ r = 0.5i(Vi7°i/' - V;7V) (2.15) 

and 

p = 2XKx + YKy - K, p^K. (2.16) 
Note that the equations of g-essence (2.9)-(2.14) can be rewritten as 



- p 


= 0, 


(2.17) 


2H + 3H^+p 


= 0, 


(2.18) 


{a^Kx^)t - a^'K^ 


= 0, 


(2.19) 


{a'KYi;^)t-2iK^{j''ij), 


= 0, 


(2.20) 


{a''KYi;f)t + 2iK4{pf), 


- 0, 


(2.21) 


p + 3H{p + p) 


= 0. 


(2.22) 



Also we present the useful formula 

KyY = 0.biKy{iPj°iP - 7/^7V) = -0.5iK^,^j + K^tJ)) (2.23) 
and the equation for u ~ 

[In {ua'^Ky^u = -iKy\TpfK^ - K^j°^P). (2.24) 

Finally, we note that some exact solutions of g-essence (2.9)-(2.14) were presented in [I7]-[I1] (see, 
also [11]- [H] in which some integrable aspects of the FRW models were considered). 

2.1 Purely kinetic g-essence 

Let us consider the purely kinetic g-essence, that is, when K = K{X, Y). In this case, the system 
(2.9)-(2.14) becomes 



-p 


= 0, 


(2.25) 


2H + 3H^ +p 


= 0, 


(2.26) 


a^Kx4> - cr 


= 0, 


(2.27) 


a^'Kyij:] - 


= 0, 


(2.28) 


a^Kyr,^ - <;* 


= 0, 


(2.29) 


p + 3H{p + p) 


= 0, 


(2.30) 
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where a (<;) is the real (complex) constant. Hence we immediately get the solutions of the Klein- 
Gordon and Dirac equations, respectively, as 



Also the following useful formula holds: 

X = ^^ or Kx^^^. (2.32) 

2.2 K-essence 

Let us now we consider the following particular case of g-essence (2.1): 

K = Ki= Ki{X,(j)) (2.33) 

that corresponds to k-essence. Then the system (2.9)-(2.14) takes the form of the equations of 
k-essence (see e.g. [3]-[B]) 

m'^-pk = 0, (2.34) 

2H + 3H^+Pk = 0, (2.35) 

Kix4>+iKix+3HKix)<p-Ki^ = 0, (2.36) 

Pk+3H{pk+Pk) = 0, (2.37) 

where the energy density and the pressure are given by 

Pk^2XKix-Ki, Pk^Ki. (2.38) 

As is well known, the energy-momentum tensor for the k-essence field has the form 

= KxV f^(jNu(j) - gp^ijK = 2KxXui^_,uiu ~ Kg^^ = {pk + Pk)ui^ui^ -pfcg^^. (2.39) 

It is interesting to note that in the case of the FRW metric (2.7), purely kinetic k-essence and 
F(T)-gravity (modified teleparallel gravity) are equivalent to each other, if a = e ^ IH]~[5]- 

2.3 F-essence 

Now we consider the following reduction of g-essence (2.1): 

K^K2 = K2{Y,4,,^) (2.40) 
that is, f-essence |13j . The energy- momentum tensor for the f-essence field has the form 

2 S S 

- gi_,yK = KYYu2t^U2u - Kg^y = {pf + Pf)u2^U2v - Pfgi^iy (2-41) 
For the FRW metric (2.7), the equations of the f-essence become [13] 

3H^~Pf = 0, (2.42) 

2H + 3H^+Pf = 0, (2.43) 

K2Yi' + 0-H^HK2Y + K2Y)yJ~il''K2^ = 0, (2.44) 

K2Y4' + 0.5{3HK2Y + K2Y)^ + iK2^l'' = 0, (2.45) 

Pf + 3Hipf+pf) = 0, (2.46) 

where 

Pf=YK2Y-K2, Pf = K2. (2.47) 
Some properties of f-essence were studied in [52]-|53). 
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3 Model 

In this section, we consider the action (2.1) with the fohowing particular g-essence Lagrangian: 

K = aiX + aaX" + a3^i((/)) + ^lY + /SaF™ + /33V'2(V^, V-) + vUii<P)U2i^, (3.1) 

where aj,l3j,ri are some real constants. As the search for exact solutions of the coupled system 
of differential equations (2.9)-(2.14) for the g-essence Lagrangian (3.1) is a very hard job, let us 
simplify the problem e.g. as V2 — V2(u),/32 — 0,Ui — (j),U2 — u. Then the system (2.9)-(2.14) 
takes the form 

ZH^ - p = 0, (3.2) 
2H + 3H^+p = 0, (3.3) 



[3F+(ln(ai+a2nX"-^))t]</)--:p^^— ^ - 0, (3.4) 

i}j + l.5Hi,-i/3-'j''{l33V2ui' + vH^) - 0, (3.5) 

ijj + 1.5Hi; + ip-\p3V2ui' + V^i'h'' - 0, (3.6) 

p + 3Hip + p) = 0, (3.7) 



where 



p = aiX + a2i2n-l)X'" -asVi- /33V2-v<l)u, (3.8) 
p = aiX + a2X" + asVi + PiY + P3V2 + Vfl^u. (3.9) 

Hence and from (2.23)-(2.24) we get 

p + p = 2aiX + 2a2nX'' + PiY ^ ~2H, (3.10) 
PiY = -iP3V2+V(^)u, (3.11) 

(3.12) 



4 Solution 

In this section, we want to present the exact solution of the system (3.2)-(3.7). But first note that 
for six unknown functions a, </), "0, i/;, Vi, V2 we have five differential equations (3.2)-(3.6) so that 
we need one more equation (see e.g. [47]-[49]). Such an equation wc take to be 

a = ao0^ (4.1) 

where oq — cq^q k = X/5. Then we obtain the following solution 

a = aot^, (4.2) 
= 00^', (4.3) 

"0 ^ 

= ^Skrx^'''^ (^ = 3,4), (4.5) 
"0 



where Cj obey the condition 
and 



ciP + |c2p-|c3p-|c4p (4.6) 



^ f ^2(5-l)+3A+l , 

Awo[2(5-l) + 3A + l] 



uo[2n((5- 1) + 3A + 1] Mo(3A-l) 
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This solution is correct if 

d = 2-3X (4.8) 

or 

The corresponding potentials take the form 

2(a-l) 2n(i5-l) 5-3X 

2(1-^) 2.i(l-iS) _ 3X-S 

u \ '■'^ / u \ ( u 



V2{u) = — + 92 — + 93 — + 94 — - 03/33-^^10(4.11) 

V"o/ V"o/ V"o/ \"o/ 

Here Vio = const, uq = ccq ^/Jf^, 

, c^i<?^0§(<?-l + 3A) 

= 2a,{S-l) ' ^^-^'^ 

, a2<5^"0g"[(2n~l)(^-l) + 3A] 

- 2"a3(<5-l) ' ^'-^'^ 

^3 = -^f^ (4.14) 
a3(d - 3A) 

and 

qi = O.5ai<5202_a3/^^ (4.15) 

q2 = a2{2n-l)2~"S^"cj)l" -ash, (4.16) 
Qa = -"3^3 - V(l^ouo, (4.17) 
Qi = -3A2. (4.18) 



5 Unified scalar— fermionic Chaplygin gas model of dark en- 
ergy and dark matter from g-essence with Yukawa inter- 
actions 

The most popular models of dark energy and dark matter such as e.g. the ACDM and a quintesse- 
CDM model assume that dark energy and dark matter are distinct entities. Another interpretation 
of the observational data is that dark energy and dark matter are different manifestations of a 
common structure. The first definite model of this type was proposed in [5T], based upon the 
Chaplygin gas, an exotic perfect fluid obeying the equation of state (EoS) 

P---, (5.1) 
P 

which has been extensively studied for its mathematical properties [5D]. The general class of 
models, in which a unification of dark energy and dark matter is achieved through a single entity, 
is often referred to as quartessence. Among other scenarios of unification that have recently been 
suggested, interesting attempts are based on k-essence. In this section we extend this scenario to 
the g-essence model with Yukawa interactions which gives us the Chaplygin gas unified model of 
dark energy and dark matter. To do it, let us consider the g-essence model given by the following 
scalar-fermionic DBI Lagrangian: 

K = U^l + ViX + V2Y^, (5.2) 

where in general U — C/(0, V',V'), Vi = Vi{(j),'ip,^/j), V2 — ^2(0, '0,V'). Note that the g-essence 
model (5.2) is constrained in two particular cases: i) the scalar DBI model as U — U{(t)), V\ = 
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Vi{(t>), V2 = 0; ii) the fermionic DBI model as U = U{ip,'ilj), Vi =0, ¥2^ V2(i^,ip). Substi- 
tuting the expression (5.2) into (2.16) we get 



where we assume that 



p = U^/l + ViX + V2Y^, 
U 



(5.3) 
(5.4) 

(5.5) 



It is the g-essence model with Yukawa interactions U = rjcfiu. These equations give the following 
EoS: 

p = - — . (5.6) 

P 

It is the Chaplygin gas model |51' but with the variable function U (Yukawa interactions). 
From (5.6) and (2.14), we get 



p = a 



-3 



; j U^{a)a 



0.5 



,-0.5 



C+ / U'^{z)dz 



0.5 



(5.7) 



where C = const, 
parameter: 



a^. From these formulas we obtain the following expression for the EoS 



_ zd\np^ _ d\np^ _ zU^ _ zh 

^ ~ ~ ^ dlnz ^ ^C + JU^{z)dz ^ ^C + h 



^-z[\n{C + h)] , (5. 



where h = J U^{z)dz, h = dh/dz. In principle, now it is not difficult to construct solutions of 
the g-essence equations corresponding to the different expressions for U. Here we just present the 
expressions for the energy density and pressure. Consider some examples. 

i) Let U — p, — const. Then from (5.7) and (5.6) we obtain the expressions for the energy 
density and the pressure 



-0.5 



;a + p^ 



0.5 



and 



c 



-z 



0.5 



(Ca-6 + ^2)0.5' 



respectively. The EoS parameter is given by 



p^z^''+^ 



C 



-z 



(5.9) 
(5.10) 

(5.11) 



So this example corresponds to the usual Chaplygin gas [HI]. As is well known, in this case, for 
small a (a^ << Cp~^'^), the energy density and the pressure take the forms, approximately, 

/?«C"'^a-^ pwO (5.12) 

with cj = 0, which corresponds to a matter-dominated universe. For a large value a, it follows that 

p~p, p^-p (5.13) 

that is, id — —1, which corresponds to a dark energy-dominated universe. So this simple and 
elegant model smoothly interpolates between a dust dominated phase, where p « C°'^a~^, and 
a de Sitter phase, where p « — p, through an intermediate regime described by the EoS for stiff 
matter, p = p. 
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ii) Now let us consider the case when U = jiz^ , where fx and v are some real constants. In this 
case, Eq.(5.7) gives the following expression for the energy density: 



-0.5 



C- 



2v + l 



0.5 



Ca 



n 0.5 



12i/ 



2zy+ 1 



(5.14) 



Here V must be negative, because otherwise, a — > oo implies p — > oo, which is not the case for 
expanding Universe. The equation (5.6) gives the expression for the pressure. 



The EoS parameter is given by 



0.5 



0.5 • 



(5.15) 



12l/ 



2v+l ' 2v+l 

The deceleration parameter q has the expression 



Ca-6(i+2.) + ^ 



a 



2M^(3t/+2) 121. 
2I/+1 



2/5 2(Ca-6 + 2^ai2-) 
For accelerating universe, q must be negative i.e., a > 0. Hence we have 

^-6(2.+i) < 



2/x2(3;^ + 2) 



C(2i/ + 1) 



(5.16) 



(5.17) 



(5.18) 



This means that for a small value of the scale factor we have a decelerating universe while for large 
values of scale factor we have an accelerating universe. The transition between these two phases 
occurs when the scale factor is equal to 



a = Oc = — 



C{2v + 1) 



-6(2ix + l) 



(5.19) 



6 Unified scalar— fermionic generalized Chaplygin gas model 
of dark energy and dark matter from g-essence with 
Yukawa interactions 

In the previous section we constructed the unified scalar-fermionic Chaplygin gas model of dark 
energy and dark matter using the particular g-essence with Yukawa interactions having the scalar- 
fermionic DBI Lagrangian form (5.2). In this section we extend these results of the previous 
section for the scalar-fermionic generalized Chaplygin gas case. For this purpose, we consider the 
particular g-essence model with Yukawa interactions which has the following scalar-fermionic DBI 
Lagrangian form 

K=U{l + ViX^ +V2Y^Y, (6.1) 



where 

U = r](j>u, = V^{4>, V), V2 = V2{4>, V)- 
Substituting this expression into (2.16) we get 

p = -U{l + ViX^ +V2Y^Y-^. 



These equations give 



i-U) 



a+1 



A 



(6.2) 

(6.3) 
(6.4) 

(6.5) 
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where a — n(\~n) ^, A — (— C/)"+^. It is the scalar-fermionic generalized Chaplygin gas model 
[5T| . From (2.16) and (6.5) we get 



3(1 + a) / (-C/)"+ia3(i+")-ida + c 



(6.6) 



where C = const, 
expression 



= a'^^^^'^K At the same time, for the pressure we obtain the following 



P = -{-U) 



a+l 



Zl + a 



C+ / (-[/)"+Mz 



These formulas give the following expression for the EoS parameter: 

z(-f7)"+i 



C + J{-U)°'+^dz' 



(6.7) 



(6.8) 



Now we can consider different types solutions of the g-essence equations. For example, if for 
simplicity, we consider the case a; = const then from (6.8) we obtain the following expression for 
the function U : 

U = -A^a-3(i+^), (6.9) 

where Aq is an integration constant. Let us consider some examples: i) if w = — 1 (the de Sitter 

1 1 
case) then U = — Aq+° ; ii) if = (the dust case) then U — — Ag+°a~^; iii) ii oj — 1 (the stiff 

1 

matter case) then U — — Aq+" and so on. Now the construction of solutions of the Friedmann, 
Klein-Gordon and Dirac equations is a formal problem so we omit it here. 



7 Scalar— fermionic DBI generalization of the g-essence 
model with Yukawa interaction (3.1) 

Our aim in this section is to construct the scalar-fermionic DBI generalization of the g-essence 
model with Yukawa interaction (3.1). To do it, we note that the model (3.1) is some approximation 
of the following scalar-fermionic DBI model: 



K 



eA + A-^[2aiX + 2a2X" + 2/3iF + 2/32^™] -l}+V, (7.1) 
where e = +1, A = const and 

V = asVii^b) + l33V2{^,i;) + TiUi{^)U2{i>,^). (7.2) 



From (7.1) we get 

P ^ 



{ v/1 + A-^ [2aiX + 2a2X" + 2/3iy + 2/32 y™] - l} + V", (7.3) 

^ v/1 + A'^[2aiX + 2a2X» + 2^^ + 2/^2 F"] 

The study the system of equations of g-essence (2.9)-(2.14) with the expressions for the pressure 
and the energy density given by (7.3)-(7.4) is a very hard job. Let us simplify the problem. Let 
0^2 = /?! = and m — 2. Then equations (7.3)-(7.4) take the form 

p - eA { + A-^2aiX + 2/32^2] - l} + (7.5) 

p = +eA-V (7.6) 

v/l + A-i[2aiX + 2/32^2] 

which corresponds to the EoS 

P = ^—--eA + V. (7.7) 

p — eA+ V 
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The corresponding EoS parameter is given by 



p{p-€A + V) p eA-V\p p-eA + vJ' 

Now let us consider two particular limit cases when p >> eA ~ V or p << eA — V . 

i) First we consider the case when p » eA — V . Let us rewrite the equation (5.7) as 

Hence, in the first approximation, we get 

42 

p = eA + V (7.10) 

P 

which corresponds to the Chaplygin gas case 151'. Now let us consider the second approximation. 
Then from (7.9) we obtain 

A^ ( tA-V\ , A^ A^eA-V) 

p = — -« 1 + ]-eA+V = ^-^ ^--eA + V. 7.11 

"^p-eA + Vpy p J p p2 



ii) Now let us consider the second case when p << cA — V. Then from (7.7) we get 

12 



where 



eA-V 

In this case, the EoS parameter is given by 



^9 = Tir^-^^ + y- (7-13) 



(eA-V)^ p 

8 Conclusion 



(7.14) 



Recently, the so-called g-essence has been proposed to be a candidate of dark energy. It has the 
non-standard kinetic scalar and fermionic terms and has interesting properties. We expect that 
g-essence has rich properties and when it is used in cosmology, its rich properties could have some 
interesting consequences. In this work, we found the exact power law solution of the particular 
g-essence model with the Yukawa scalar-fermionic interactions. The corresponding scalar and 
fermionic potentials are presented. Of course, there might be other exact solutions of the g-essence 
different from the power law solution. Anyway, our results obtained in the present work showed 
that the g-essence with the Yukawa interactions can describes the accelerated expansions of the 
universe. In fact, we have shown that some particular g-essence models with Yukawa interactions 
correspond to the usual and generalized Chaplygin gas unified models of dark energy and dark 
matter. Also we presented some scalar-fermionic DBI models which again can describe the unified 
dark energy-dark matter system. 
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